rm(list = ls(all = TRUE)) #clear workspace

library(foreign)
library(survey)
options(survey.lonely.psu = "adjust")

setwd("C:/Users/martin.romero/Dropbox/ENVIPE/2010")
#setwd("~/Dropbox/ENVIPE/2011/Base_de_datos")

################################################
#Pregunta al entrevistado
################################################

data <- read.dbf("TDELITOS2.dbf")

#Ponderador y diseno
data$FAC_ELE <- as.numeric(as.character(data$FAC_ELE))
design.per <- svydesign(id=~UPM,strata=~EST, weights=~FAC_ELE, data=data)

#4.1 Mencion� que en 2009 usted fue v�ctima de alg�n delito. 
#�Podr�a decirme de cu�l o cu�les delitos se trat�?
#Secuestro o secuestro expr�s
#1 S�
#2	No
#b	blanco

table(data$DP4_1_7, useNA="ifany") #41 casos
data$DP4_1_7 <- factor(data$DP4_1_7, levels=c(0,1))
data$DP4_1_7[is.na(data$DP4_1_7)] <- 0
table(data$DP4_1_7, useNA="ifany") # check
design.per <- svydesign(id=~UPM,strata=~EST, weights=~FAC_ELE, data=data) # ...again
foo <- svytotal(~DP4_1_7, design=design.per, na.rm=TRUE) #67,677
print(foo)
confint(foo) #35,036-100,317

#4.2 �Cu�ntas veces le sucedi�?   					
table(as.numeric(as.character(data$DP4_2_7)))
sum(as.numeric(as.character(data$DP4_2_7)), na.rm=TRUE) #43
foo <-svytotal(~as.numeric(as.character(DP4_2_7)), design=design.per, na.rm=TRUE) #70,256
print(foo)
confint(foo) #36,640-103,871
cv(foo)
